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Abstract 

The Eulerian system of dynamic equations for the ideal (nondissipative) 
fluid is closed but incomplete. The complete system of dynamic equations 
arises after appending Lin constraints which describe motion of fluid particles 
in a given velocity field. The complete system of dynamic equations fo the 
ideal fluid can be integrated. Description in terms of hydrodynamic potentials 
(DTHP) arises as a result of this integration. The integrated system contains 
indefinite functions of three arguments, which can be expressed via initial and 
boundary conditions. The remaining initial and boundary conditions for the 
integrated system can be made universal (i.e. the same for all fluid flows), 
and the resulting system of equations contains full information about the 
fluid flow including initial and boundary conditions for the fluid flow. Some 
hydrodynamic potentials appear to be frozen into the fluid, and the Kelvin's 
theorem on the velocity circulation can be formulated in a contour-free form. 
Description in terms of the wave function (DTWF) appears to be a kind of 
DTHP. Calculation of slightly rotational flows can be carried out on the basis 
of DTHP, or DTWF. Such a description of a rotational flow appears to be 
effective. 



1 Introduction 

Dynamic equations for an ideal (nondissipative) fluid are written conventionally in 
the form: 

f + V(pv)=0 (1.1) 
^ + (vV)v = --Vp. p = p=^ (1.2) 

— + (vV)5 = (1.3) 

where dependent variables p and v = {v^,v'^,v^}, S are respectively the fluid mass 
density, the fluid velocity and entropy per unit mass considered as functions of 
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independent Eulerian variables x — {t, x}. p is athc pressure, and E — E{p, S) is an 
internal energy of an unite mass considered as a function of p and S. The internal 
energy E = E{p, S) is an unique characteristic of the ideal fluid. 

The system of hydrodynamic equations (1.1)-(1.3) is a closed system of differ- 
ential equations which has an unique solution inside some space-time region fl, pro- 
vided dependent dynamic variables p and v = {v^,v'^,v^}, S are given as functions 
of three arguments on the space-time boundary F of the region Q. Being closed, the 
system (1.1)-(1.3) describes nevertheless only momentum-energetic characteristics 
of the fluid. Motion of the fluid particles along trajectories is described by so called 
Lin (1963) constraints 

§ + (vV)C = 0, (1.4) 

where quantities ^ = ^(t, x) = {^^(^jx)}, a = 1,2,3 label fluid particles. They 
will be referred to as particle labeling (curvilinear Lagrangian coordinates). If the 
equations (1.4) are solved and ^ is determined as a function of (t,x), the finite 
relations 

^(i, x) = Cin = const (1.5) 

describe implicitly a fluid particle trajectory and a motion along it. 

The system of eight equations (1.1) -(1.4) forms a complete system of dynamic 
equations describing a fluid, whereas the system of five equations (1.1)-(1.3) forms 
a curtailed system of dynamic equations. The last system is closed, but to be a 
complete system, it must be supplemented by the kinematic equations 

dx 

— = v(i,x), x = x(t,0 (1.6) 

where v(t, x) is a solution of the system (1.1)-(1.3). Three equations (1.6) are 
equivalent to (1.4), because any solution ^ = ^{t, x) of (1.4) is a set of three integrals 
of equations (1.6). 

There is a lack of understanding of the fact that the Euler system (1.1)-(1.3) is a 
curtailed one and that the equations (1.4) describe a motion of fluid particles. Even 
Lin (1963) who introduced equations (1.4) considered them "as the condition for the 
conservation of the identity of particles" , but not as kinematic equations describing 
a motion of fluid particles in a given velocity fleld. Betherton (1970) investigated the 
Lin constraints and showed a connection between them and the Kelvin's theorem on 
the velocity circulation. But the very important statement that the Lin constraints 
(1.4) are a kind of kinematic equations (1.6) was not mentioned directly. In general, 
we failed to flnd in literature any reference to equations (1.4) as necessary kinematic 
equations, although this fact seems to be evident. Apparently, it is connected with 
the common belief that the Lin constraints (1.4) are useless, if one is interested only 
in the velocity fleld of the flow. 

We should like to show a necessity of these equations in a very simple example, 
modelling a situation in the hydrodynamics in a grotesque form. Let us consider a 
particle moving in uniform gravitational field g =const. Dynamic equations have 
the form 

V = g, X = V, g = const (1.7) 
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where x and v are functions of t describing respectively position and velocity of 
the particle. First three equations (1.7) constitute a closed subsystem [analog of 
the Euler system (1.1)-(1.3)] of the full system of six equations. Equations of the 
subsystem can be solved independently of the remaining equations. But it does 
not mean that this closed subsystem may be considered as a system of dynamic 
equations describing a particle, even if we are interested only in momentum-energetic 
properties of the particle. Such important characteristic of the particle as the energy 
integral — gx = const cannot be derived on the base of only first three equations 
(1.7). This integral is a sum of integrals 

{Vaf/2 - QaXa = C„ = COUSt, Va = , = 1, 2, 3 (1.8) 

These integrals are analogs of integrals of the complete system (1.1)- (1.4) [see below 
(1.11)]. This example shows that a closed system of equations and a complete system 
of dynamic equations is not the same. The Lagrangian formulation of hydrodynamic 
equation includes equations (1.4) in the form (1.6) automatically. Lagrangian for- 
mulation is equivalent to the Eulerian formulation (1.1)-(1.3), provided equations 
(1.4) are appended to it. 

Equations (1.1), (1.3) can be integrated on the basis of (1.4) in the form 

5(t,x) = 5o(0 (1-9) 

p(t,x)=po(l) g(^,^^,^^3) -Po(0^ (1-10) 

where Sq{$,) and po(^) ^'I'c arbitrary integration functions of the argument ^. These 
functions can be determined from the initial (and boundary) conditions. Three 
equations (1.2) also can be integrated on the basis of (1.4). These integrals have the 
form 

v(i, x) = u{<p, ^, 77, S) = V<p + ^"(^) VCa - V^S, (1.11) 

where g{^) — {g'^i^)}, ol — 1, 2, 3 are arbitrary integration functions of argument ^, 
and </?, 77 are new dependent variables satisfying dynamic equations 

^ + u((^, 77, S)V^ - \[n{^, 77, S)f + ^ = (1.12) 

^ + u((^,^,77,5)V77=-||. (1.13) 

If five dependent variables ^, 77 satisfy the system of equations (1.4), (1.12), (1.13), 
the five dynamic variables 5, p, v (1.9)-(1.11) satisfy dynamic equations (1.1)-(1.3). 
Indefinite functions (?"(^) can be determined from initial and boundary conditions 
in such a way that the initial and boundary conditions for variables ^, r] were 
universal in the sense that they do not depend on the fluid flow. 

According to (1.10), (1.11) the physical quantities p, v are obtained as a result of 
differentiation of the variables </?, ^, and the variables </?, ^, 77 can be regarded as 
hydrodynamic potentials. These potentials appear in the Hamilton fluid dynamics 
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(Salmon, 1988) as dependent variables. The hydrodynamic potentials arise by a 
natural way. They associate with the name of Clebsch (1857, 1859) who introduced 
these quantities for the incompressible fluid. Such quantities as g°'{$) also appear in 
the Hamilton fluid mechanics (Salmon, 1988), but they appear as dependent variables 
(Lagrage invariants) satisfying dynamic equations of the type (1.4). They also are 
regarded as hydrodynamic potentials. Note that in the Hamilton fluid mechanics 
the quantities g°' are considered simply as dependent variables, but not as indefinite 
functions of |, arising as a result of integration, although corresponding dynamic 
equations for g"" can be integrated easily. 

It should distinguish between the integration and a change of variables which 
does not contain arbitrary functions explicitly. For instance, let us substitute g"{^), 
q; = 1, 2, 3 by new dependent variables A"', a — 1,2,3, imposing on them constraints 

-^ + u{<p,^,ri,S)VA''^0, a = 1,2, 3. (1.14) 

Then one has instead of (1.11) 

V{t, X) = U((^, ^, T), S) = V^ + ^l^VCa - V^S, (1.15) 

The eighth order system (1.4), (1.12)-(1.15) arises instead of the fifth order 
system (1.4), (1.11)-(1.13). This system is not an integrated one, because it has 
the higher order and contains only two arbitrary functions (1.9), (1.10). Different 
modifications of such kind systems were derived (Salmon, 1988). They cannot be 
considered as integrated systems. Of course, one can easily integrate equtions (1.14) 
on the basis of (1.4), and return to (1.11). But until this is not made and the 
number of dynamic variables is not decreased, the system cannot be considered as 
integrated. 

The integration of the complete system (1.1)-(1.4) and some corollaries of this 
integration correlates with the Hamilton properties of the ideal fiuid [Harivel, (1955); 
Eckart, (1960); Seliger and Whithem, (1967), Salmon, (1988); Zakharov and Kuznetsov 
(1997)]. It is connected with the fact that the curtailed system (1.1)-(1.3) is not a 
Hamiltonian system in itself, whereas the complete system (1.1)-(1.4) is a Hamil- 
tonian one that can be easily seen in the example (1.7). Constructing Hamiltonian 
mechanics of the ideal fluid, one uses (implicitly or explicitly) the Lin constraints 
(or part of them). It is this expansion of the curtailed system (but not Hamiltonian 
properties) that is important for integration and derivation of other useful results. 
To show this, the Hamiltonian technique and Hamiltonian properties of the ideal 
fluid will not be used at all. 

Note that the curtailed system (1.1)-(1.3) has the same order as the integrated 
system (1.4), (1.12), (1.13), but it takes into account neither initial conditions, 
nor kinematic equations (1.6). The fact that the ideal fluid considered as a dynamic 
system admits both the curtailed system (1.1)-(1.3) and the integrated system (1.4), 
(1.12), (1.13) is connected closely with the group of the relabeling transformation 
(relabeling group) 

^a^L^U^), = det II 0, a, (3 ^1,2,3 (1.16) 
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<^ = ^0 ^ lo = <^ = lo(6) + aoii), d^o/d^o > (1.17) 
where ^ = {^o, ^} are curvilinear Lagrangian coordinates in the space-time, ^ = 
{Co, i} is another system of curvihnear Lagrangian coordinates, and qq are arbi- 
trary functions of ^. is arbitrary function of ^o- Co is a temporal coordinate, and 
^ are spatial ones. 

The relabeling group is a symmetry group of the fluid considered as a dynamic 
system. This circumstance admits to integrate the complete system (1.1)-(1.4). Any 
special particle labeling is unessential from physical viewpoint. It is a reason why 
several equations (1.1)-(1.3) of the complete system form a closed system describing 
conservation laws. The relabeling group is used in hydrodynamics comparatively 
recently Eckart (1938, 1960), Calkin (1963), Bretherton (1970), Friedman and Schutz 
(1978), Salmon (1982), Zacharov and Kuznetsov (1997) and others. 

The integrated system (1.4), (1.12), (1.13) looks more comphcated, than the 
curtailed system (1.1)-(1.3). It is quite natural, because the integrated system is a 
complete system which contains five indefinite functions po, Sq, g°' describing initial 
and boundary conditions. 

There is a common belief that the curtailed Eulerian system (1.1)-(1.3) is suffi- 
cient for calculating fluid flows, and a use of the complete system (1.1)-(1.4) is not 
necessary in most of cases. In general, there is a puzzling question of such a kind. 
If both the Euler system (1.1)-(1.3) and the integrated system (1.4), (1.12), (1.13) 
have the same order and the same number of dependent variables, why should one 
consider the integrated system (1.4), (1.12), (1.3) which looks more complicated, 
than the Euler system? The answer is as follows. Although the Lin constraints de- 
scribe mainly a motion of fluid particles in the given velocity field, nevertheless the 
integrated system contains additional information which is necessary for a calcula- 
tion of rotational fluid flows. This additional information concerns fluid properties 
described by the Kelvin's theorem on the velocity circulation. This theorem is an 
attribute of the complete system (1.1)- (1.4), because it refers to the contour con- 
nected rigidly with the fluid particles moving according to dynamic equations (1.6), 
or (1.4). In general, the Kelvin's theorem cannot be formulated only in terms of 
the velocity fleld. On the other hand, it is clear that the fluid properties (vorticity 
frozen into fluid) described by the Kelvin's theorem are very important for calculat- 
ing rotational flows, whereas for irrotational fiows the constraints imposed by the 
Kelvin's theorem degenerate into identities which are fulfilled automatically. 

The conventional theory of fluid flows is based mainly on the Euler system (1.1)- 
(1.3) which does not take into account constrains of the Kelvin's theorem. Theory 
of irrotational flows has been developed well enough, whereas that of rotational 
flows has been developed much shghtly. Apparently, it is connected with the fact 
that the Euler system (1.1)- (1.3) does not contain enough information on the fluid 
properties, and proper calculations of rotational flows are impossible in the scope 
of the curtailed Euler system (1.1)-(1.3). Especially it concerns strongly rotational 
(turbulent) flows. It will be shown in this paper that slightly rotational flows of 
incompressible fluid can be calculated on the basis of the integrated system (1.4), 
(1.12), (1.13) (DTHP or DTWF), where integrals (1.11) are taken into account. 

The integrals (1.11) can be interpreted as a contour-free form of the Kelvin's 
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theorem in application to the perfect fluid. Indeed, let us multiply (1.11) by dx and 
integrate along a closed path C One obtains 



Cc = / vdx ^jg-{i)di^ + ^ri^^d^a (1.18) 

where is a mapping of the contour C in the x-space onto the ^-space of labels 
This mapping 

C^C^. ^^x(t,^), x^^(i,x) (1.19) 

depends on time t. Let the contour C be coupled rigidly with fluid particles and 
move with the fluid. It means that is flxed, x(t,^) satisfles (1.6), and the shape 
of C depends on time according to (1.19), (1.6). Let the flow be homoentropic 
{S{^) =const identically), or the contour jC^ lies on the surface S{^) =const., the 
second integral in rhs of (1.18) vanishes, and the circulation Cl does not depend 
on time. 

Let now the inflnitesimal contour £ be a parallelogram made up by two inflnites- 
imal vectors (ixi, (ix2 and be made up by two inflnitesimal vectors d^i, d^2- Then 
using Stokes's theorem, one derives from (1.18) 

u^dS = n{i)dS^, w = V X V, n{i) ^ {Oa(C)}, a = 1, 2, 3 (1.20) 

= « = 1, 2, 3 dS = dxi X dx2, dS^ = d^i X d^2 

^ap _ _ dg^ a 3-123 (121) 

where Sap-y is the Levi-Chivita pseudotensor, and a summation is produced (1-3) 
over repeated Greek indices. dS and dS^ are infinitesimal area of contours C and 

respectively. The relation (1.20) is a local (or contour-free) form of the Kelvin's 
theorem. The scalar wdS conserves and does not depend on time, although a; and dS 
individually depend on time. For any irrotational flow the relation (1.20) degenerates 
into identity, and Kelvin's theorem may be ignored. 

Note that f2(^) may be regarded as a "frozen vorticity", because it depends on 
time t only via ^, and f2 = provided ^ = x. r2(^) is a scalar in the x-space, and 
it is a vector in the |-space of labels |. 

For derivation of integrated system (1.4), (1.12), (1.13) one uses a speciflc mathe- 
matical technique based on Jacobian properties. This technique permits to integrate 
dynamic equations without a use of a change of variables. In general, all results can 
be obtained, carrying out a proper change of variables in the action functional by 
means of the Lagrange multipliers. Such changes of dependent variables are pro- 
duced in the Hamilton fluid dynamics (see, for instance, Salmon, 1988). Dynamic 
equations of the type of (1.3) appear as a result of such changes of variables. These 
equations can be integrated easily on the basis of equations (1.4). Unfortunately, 
such changes of variables lead to different sets of dependent variables whose physi- 
cal meaning is unclear. In other words, using Lagrange multipliers, for a change of 
variables, one obscures logical connection between different variables. It is rather 
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difficult to understand that, integrating some equations of the type of (1.3), one 
integrates in reahty the equations (1.2). To simphfy the logical connection between 
different dependent variables and to clear their physical meaning, we prefer to inte- 
grate dynamic equations directly by means of "Jacobian technique" (sec. 2). Use 
of Jacobians in hydrodynamics has had a long history, dating back to the time of 
Clebsch (1857, 1859). It was the use of Jacobians that allowed to introduce the 
Clebsch potentials and integrate hydrodynamic equations. 

The Jacobian technique was used by many authors (Hcrivel (1955), Eckart(1960), 
Berdichevski (1983), Salmon (1988), Zacharov and Kuznetsov (1997) and many 
others). We use space-time symmetric version of the Jacobian technique which 
appears to be simple and effective. It seems that the progress in the integration of 
hydrodynamic equations is connected mainly with the developed Jacobian technique. 



2 Jacobian technique 

Let us consider such a space-time symmetric mathematical object as the Jacobian 
J ^ 51°' ^2^1 ^ det II II, = 46 ^% ^,k = 0, 1, 2, 3 (2.1) 

Here ^ = {^o,^} = {^0,^1,^2,^3} four scalar considered as functions oi x — 
{a;°,x}, ^ = ^{x). The functions {^0,^1,^2,^3} are supposed to be independent in 
the sense that J 7^ 0. It is useful to consider the Jacobian J as 4-linear function of 
variables C,i,k = dk^i, i,k = 0, 1,2,3. Then one can introduce derivatives of J with 
respect to 6,^. The derivative dJ/d^i^k appears as a result of a substitution of ^ by 
x'' in the relation (2.1). 

dJ _^{Co,■■■^^-l,x^^.^^,...i^l^ i,A; = 0,l,2,3 (2.2) 



d^i,k d{x'^,x^,x'^,x^) 
For instance 

dJ _ 9(,r'. (1.(2.^3) 



dCQ,i d{x^,x\x'^,x^)' 
This rule is valid for higher derivatives of J also. 



i = 0,l,2,3 (2.3) 



dCi,kdCs,i d{x^,x^,x'^,x^) 

(9(x^x') (9(^0,6,6,6) _ ji dx'^dx^ dx''dx\ _ni9Q (o a\ 

5(6,6)5(Aa:^a;^a;3)-^^96c?6 ^6^6^' ^' ^ " ^' ^' ^ ^'■^) 
It follows from (2.1), (2.2) that 

dx'' _ a(6,...6-l-T^6+l,••■6) _ 5(6, •••6-1, 6+1, •••6) 



^6 5(6,6,6,6) d{x^,x^,x^,x^ 

d{x'^,x\x^,x^) _ 1 dJ 



X 



5(6,6,6,6) ■^56,fe' 



i, A; = 0,1, 2, 3 (2.5) 
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and (2.4) may be written in the form 

9V I , dJ dJ dJ dJ 



i,k, I, s ^0,1,2,3 (2.6) 



The derivative dJ/d^i^k is a cofactor to the element ^i^k of the determinant (2.1). 
Then one has the following identities 

dJ dJ 
Ci,k^ = StJ, ^k,i^ = StJ, 5 = 0,1,2, 3 (2.7) 

dJ d'^ J 

^fe^ = ^fe^'^^ = Q' ^ = 0,l,2,3. (2.8) 

Here and further a summation on two repeated indices is produced (0-3) for Latin 
indices and (1-3) for the Greek ones. The identity (2.8) can be considered as a 
corollary of the identity (2.6) and a symmetry of dkdi^s with respect to permutation 
of indices k, I. Convolution of (2.6) with dk, or di vanishes also. 

Relations (2.1) -(2.6) are written for four independent variables x, but they are 
valid in an evident way for arbitrary number n -|- 1 of variables x — x^, . . . x"} 

andc = {eo,^}, I = {ei,e2,...en}. 

Application of the Jacobian J to hydrodynamics is founded on the property, that 
the fluid flux 

(97 

= i = {/} = {APv}, i = 0,l,2,3 (2.9) 

constructed on the basis of the variables ^ = {^i, ^2, Cs} satisfies Lin constraints (1.4) 
and the continuity equation 

d,f = (2.10) 

identically for any choice of variables ^, as it follows from the identity (2.8) for i = 0. 
The continuity equation (2.10) is used without approximations in all hydrodynamic 
models, and the change of variables {p, pv} ^ described by (2.9) is very important. 

In particular, in the case of two-dimensional established flow of incompressible 
fluid the variables | reduce to one variable = ip, known as the stream function. 
In this case there are only two essential dependent variables x^ = x, x^ = y, and 
the relations (2.9), (2.10) reduce to relations 

-1 • 9ij _i . dip dv^ dvy 



Deflning the stream line as a line tangent to the flux j 

dx dy 
~ ~ ) 

Jx Jy 



(2.12) 



one obtains that the stream function is constant along the stream line, because 
according to two flrst equations (2.11), ip — ip{x,y) is an integral of the equation 
(2.12). 
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In the general case, when the space dimensionahty is n and x — {x^, x^, . . . x"}, 
^ = {^0) ^ — {Ci) ■ ■ ■ Cn}, the quantities ^ = {Co}, a — 1,2, . . .n are constant 
along the hne C tangent to the flux vector j — {f}, i — 0,1, . . .n 

d, -j^ = f{x), t = 0,l,...n (2.14) 

where r is a parameter along the line C which is described parametrically by the 
equation x = x{t). This statement is formulated mathematically in the form 

^ = j%^a = m^diU = 0, a = 1, 2, ... n (2.15) 

The last equality follows from the first identity (2.7) taken for s = 0, / = 1, 2, . . . n 
Interpretation of the line (2.14) tangent to the flux is different for different cases. 
If X = {x°, x^,. . . x"} contains only spatial coordinates, the line (2.14) is a line in the 
usual space. It is regarded as a stream line, and ^ can be interpreted as quantities 
which are constant along the stream line ( generalized stream function). 

If x^ is the time coordinate, the equation (2.14) describes a line in the space-time. 
This line (known as a world line of a fluid particle) determines a motion of the fluid 
particle. Variables ^ = {^1,^2, • • • Cn} which are constant along the world line are 
different, generally, for different particles, li ^a, 0( — 1,2, . . .n are independent, they 
may be used for the fluid particle labeling. 

Thus, although interpretation of the relation (2.9) considered as a change of 
dependent variables j by ^ may be different, from the mathematical viewpoint this 
transformation means a replacement of the continuity equation by some equations 
for the labeling (or generalized stream function) ^. Difference of the interpretation 
is of no importance in this context. 
Note that the expressions 

f = '-P»(()g^-"'P''(i) al^:^^% ' i = 0,1.2.3, (2.16) 

can be also considered as four-flux satisfying the continuity equation (2.10). Here 
m is a constant and Po(^) is an arbitrary function of ^. It follows from the identity 

'^P'^^ h{x^x\x\x^) ='^ d{x^x\x\x^y ^1 = 7/^0(^.6,6^,. (2.17) 

As an example of application of the Jacobian technique, let us show that (1-10) 
satisfies (1.1) in virtue of (1.4). Let us multiply (1.4) by (1.10) and introduce new 
variables j = pv = {j^,3^,3^}- One obtains three equations 

dJ 

mpo{^)^^f},o+r^p,a = 0, 13= 1,2,3. (2.18) 

Considering (2.18) as a system of three hnear equations for a = 1,2,3 and 
resolving it with respect to j", one obtains 

dJ 

r = mpo{^)^, a = 1,2,3 (2.19) 
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It is easy to verify this, substituting (2.19) into (2.18) and using (2.7). One obtains 
that j — = {p,pv} is described by the relations (2.16) which satisfy the 

continuity equation (2.10) identically. Thus, (1.1) is satisfied by (1.10) in virtue of 



In general, equivalency of the system (1.4), (1.12), (1.13) and the system (1.1)-(1.4) 
can be verified by a direct substitution of variables p, S, v, defined by the relations 
(1.9)-(1.11), into the equations (1.1)-(1.3). Using equations (1.4), (1.12), (1.13), 
one obtains identities after subsequent calculations. But such computations do not 
display a connection between the integration and the invariancy with respect to the 
relabeling group (1.16). Besides a meaning of new variables ip, rj is not clear. We 
shall use for our investigations a variational principle. Note that for a long time 
a derivation of a variational principle for hydrodynamic equaitons (1.1)-(1.3) was 
existing as a self-dependent problem (Davydov, 1949; Herivel, 1955; Eckart, 1960; 
Lin, 1963; Seliger and Whithem, 1967; Bretherton, 1970; Salmon, 1988). Existence 
of this problem was connected with a lack of understanding that the system of 
hydrodynamic equations (1.1)-(1.3) is a curtailed system, and the full system of 
dynamic equations (1.1)-(1.4) includes equations (1.4) describing a motion of the 
fluid particles in the given velocity field. The variational principle can generate only 
the complete system of dynamic variables (but not its closed subsystem). Without 
understanding this one tried to form the Lagrangian for the system (1.1)-(1.3) as a 
sum of some quantities taken with Lagrange multipliers. Ihs of dynamic equations 
(1.1)-(1.3) and some other constraints were taken as such quantities. 

Now this problem has been solved (see review by Salmon, 1988) on the basis 
of the Eulerian version of the variational principle for the Lagrangian description, 
where equations (1.4) appear automatically and cannot be ignored. In our version of 
the variational principle we follow Salmon (1988) with some modifications which un- 
derline a curtailed character of hydrodynamic equations (1.1)-(1.3), because the un- 
derstanding of the curtailed character of the system (1.1)-(1.3) removes the problem 
of derivation of the variational principle for the hydrodynamic equations (1.1)-(1.3). 

We consider the ideal fiuid as a conservative dynamic system whose dynamic 
equations can be derived from the variational principle. This dynamic system is 
a continuous set of many identical particles moving in some self-consistent (and 
external) potential force field. The action functional has the form 



where x = {x°'{t, ^)}, a = 1, 2, 3 are dependent variables considered as functions of 
time t and of labels (Lagrangian coordinates) ^ = {^1)^2)^3}- dx./dt is a derivative 
of X with respect to t taken at fixed 4- 



(1.4). 



3 Variational principle 




(3.1) 



dx' 



dt 



a = 1,2,3 



(3.2) 
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Po($) is some non-negative weight function, and y is a potential of a self-consistent 
force field which depends on ^, x and derivatives of x with respect to m =const 
is some mass of the fluid particle. The form of the potential V will be fixed later. 
Now it is important only that V does not depend on the time derivatives of x. 

Variation of the action with respect to x generates six first order dynamic equa- 
tions for six dependent variables x, v = dx/dt, considered as functions of t and of 
independent curvilinear Lagrangian coordinates ^. It is a Lagrangian representation 
of hydrodynamic equations. 

We prefer to work with Eulerian representation, when Lagrangian coordinates 
(particle labeling) ^ = {^0,^}, C = {'Ci!'C2!'C3} S'l'c considered as dependent variables, 
and Eulerian coordinates x — {a;^,x} = {t, x}, x = {x^,x'^,x^} are considered as 
independent variables. Here is a temporal Lagrangian coordinate which evolves 
along the particle trajectory in an arbitrary way. Now the ,^0 is a fictive variable, 
but after integration of equations the ,^0 stops to be fictive and turns to the variable 
ifi, appearing in the integrated system (1.4), (1.12), (1.13). 

Further mainly space-time symmetric designations will be used, that simplifies 
considerably all computations. In the Eulerian description the action functional (3.1) 
is to be represented as an integral over independent variables x = {a;°,x} = {t,x}. 
One uses the Jacobian technique for such a transformation of the action (3.1), 

Let us note that according to (2.3) the derivative (3.2) can be written in the 
form 

.«.^.J^(|L)-\ „.,2,3. (3.3) 

Then components of the 4-flux j = = {p, pv} can be written in the form 

(2.16), provided the designation (1.10) 

' = ' = '^''^^^d^o='"''^^^d{xO,x\x^x^) ^'-'^ 

is used. 

At such form of the mass density p the four-flux j = {j*}, i = 0, 1, 2, 3 satisfies 
identically the continuity equation (2.10) which takes place in virtue of identities 
(2.7), (2.8). Besides in virtue of identities (2.7), (2.8) the Lin constraints (1.4) are 
fulfilled identically 

a = 1,2, 3. (3.5) 

Components f are invariant with respect to the relabeling group (1.16), provided 
the function po(^) transforms as follows 

Po(^) Po(^) = D po(^), D = = ^ ^ ^ (3.6) 

One has 

dJ 

poU)dtd^ = po(€)7^ dtdsi = —dtdx (3.7) 

0^0.0 rn 



m 



~dt] 




(3.8) 
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and the variational problem with the action functional (3.1) is written as a varia- 
tional problem with the action functional 

M^] -[if- pE)dtdx, E^- (3.9) 

where p = j° and j = {j^, j^, j'^} arc fixed functions of ^ = {C,o, 1} and of ^a.i = diC,a, 
q; = 1, 2, 3, i = 0, 1, 2, 3, defined by the relations (2.16). E is the internal energy of 
the fluid which is supposed to be a fixed function of p and 5*0 (C) 

E^E{p,So{^)), (3.10) 

where p is defined by (3.4) and So{^) is some fixed function of ^, describing initial 
distribution of the entropy over the fluid. 

The action (3.9) is invariant with respect to subgroup Qso of the relabeling group 
(1.16). The subgroup Qso is determined in such a way that any surface Sq(^) =const 
is invariant with respect to Qso- general, the subgroup Qso is determined by two 
arbitrary functions of ^. 

The action (3.9) generates the six order system of dynamic equations, consisting 
of three second order equations for three dependent variables ^. Invariancy of the 
action (3.9) with respect to the subgroup Qso admits one to integrate the system 
of dynamic equations. The order of the system is reduced, and two arbitrary in- 
tegration functions appear. The order of the system is reduced to five (but not to 
four) , because the fictive dependent variable stops to be Active as a result of the 
integration. 

Unfortunately, the subgroup depends on the form of the function So{$,) and 
cannot be obtained in a general form. In the special case, when 5*0 (C) does not 
depend on ^, the subgroup Qso coincides with the whole relabeling group Q, and the 
order of the integrated system is reduced to four. 

In the general case it is convenient to introduce a new dependent variable 

S = Soi^). (3.11) 

According to (3.5) the variable S satisfies the dynamic equation (1.3) 

fdiS = 0. (3.12) 

In virtue of designations (2.9) and identities (2.7), (2.8) the equations (3.12), (3.5) 
arc fulfilled identically. Hence, they can be added to the action functional (3.9) as 
side constraints without a change of the variational problem. Adding (3.12) to the 
Lagrangian of the action (3.9) by means of a Lagrange multiplier rj, one obtains 

^e[^, V:S]^ l{^-pE + rjj%S}dtdx (3.13) 

where the quantities j = are determined by (2.16), and E — E{p,S). The 

action (3.13) is invariant with respect to the relabeling group Q which is determined 
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by three arbitrary functions of ^. Three arbitrary functions of ^ appear in conse- 
quence of the integration of dynamic equations. The integrated system contains five 
first order equations for dependent variables .^o, ^, V- The dependent variable S is 
substituted by arbitrary indefinite function Sq{^). 

To obtain the dynamic equations, it is convenient to introduce new dependent 
variables j\ defined by (2.16). Let us introduce the new variables f by means of 
designations (2.16) taken with the Lagrange multiphers p^, i = 0, 1, 2, 3. Then the 
action (3.13) takes the form 

AE[p,i,tP,V,S] = {^-pE-pk[f-mpoi^)^]+vj%S}dtd^ (3.14) 

It is useful to keep in mind that four designations (2.16), introducing variables p, 
j = pv via variables ^, are equivalent to three Lin constraints (1.4) together with the 
designation (3.4), as it was shown in the end of sec.2. Addition of relations (2.16) to 
the action (3.13) as side constraints is equivalent to the addition of relations (1.4), 
(3.4) considered as side constraints. 

For obtaining dynamic equations the variables p, j, rj, S are to be varied. Let 
us eliminate the variables Pi from the action (3.14). Dynamic equations arising as a 
result of a variation with respect to have the form 

—— = CaP = -m9fc Po(^) Pi\ + m pk = 0, a = 1, 2, 3 

(3.15) 

where are linear operators acting on variables p — {pi}, i — 0,1,2,3. This 
equations can be integrated in the form 

Pi = 9'{^o)diCo + 9''U)diCa, ^ = 0, 1, 2, 3, (3.16) 

where is some new variable (temporal Lagrangian coordinate), g°'{^), a = 1,2,3 
are arbitrary functions of the label ^, g^{io) is an arbitrary function of ^o- The 
relations (3.16) satisfy equations (3.15) identically. Indeed, substituting (3.16) into 
(3.15) and using identities (2.6), (2.7), one obtains 



-mdk { po(0 



+ m^^Jg'{Co)^0, a = 1,2, 3, 



(3.17) 

Differentiating braces and using identities (2.8), (2.7), one concludes that (3.17) is 
an identity. 

Setting for simplicity 

dk'P ^ g%^o)dkU A; = 0,1, 2, 3 (3.18) 

one obtains 

Pk^dkcp + g''{^)dk^a, A; = 0,1, 2, 3 (3.19) 

Substituting (3.19) in (3.14), one can eliminate variables Pi, i = 0,1,2,3 from 
the functional (3.14). The term g°'{^)dk^adJ/d^o^k vanish, the term dkcpdJ/d^o^k 
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gives no contribution into dynamic equations. The action functional takes the form 
^g[p, j, ^, rj,S]^ Ji^-pE- j^idkV + g''{i)duio. - vdkS]}dtdx (3.20) 

where g'^l^) are considered as fixed functions of ^ which are determined from initial 

conditions. Varying the action (3.20) with respect to (f, r}, S, j, p, one obtains 
dynamic equations 

Sip : dkf = 0, (3.21) 

S^^: n''f'Mip = 0, a = 1,2,3, (3.22) 
where is defined by (1.21) 

5ri : j'dkS = 0, (3.23) 

SS : fd,rj = -p—, (3.24) 
5i: v=j/p = + (?"($) V^a-r/V^, (3.25) 

~ ^"'^ " 9''i^)doL + vdoS = 0, (3.26) 

Deriving relations (3.22), (3.24), the continuity equation (3.21) was used. It is easy 
to see that (3.22) is equivalent to (1.4), provided 

det 11 Q""^ II 7^ (3.27) 

Then the equations (3.23) and (3.21) can be integrated in the form of (1.9) 
and (1.10) respectively. Equations (3.24) and (3.25) are equivalent to (1.13) and 

(1.11) . Finally, eliminating c^o^a and OqS from (3.26) by means of (3.22) and (3.23), 
one obtains the equation (1.12) and, hence, the system of dynamic equations (1.4), 

(1.12) , (1.13), where designations (1.9)-(1.11) are used. 

The curtailed system (1.1)-(1.3) can be obtained from equations (3.21)-(3.26) 
as follows. Equations (3.21), (3.23) coincide with (1.1), (1-3). For deriving (1.2) let 
us note that the vorticity u> = V x v and v x u> are obtained from (3.25) in the form 

u, = V X V = Jo"'^VC/3 X VCa -V7]xVS (3.28) 

V X w = ^^"^V^^(vV)Ca + V5(vV)77 - V77(vV)5 (3.29) 

Let us form a difference between the time derivative of (3.25) and the gradient 
of (3.26). Ehminating Jl"^9o^a, Q^S and d^i] by means of equations (3.22), (3.23), 
(3.24), one obtains 

-n"^Ve/3(vV)Ca + Vr]{vV)S - VS{vV)ri = (3.30) 
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Using (3.28), (3.29) the expression (3.30) reduces to 



-.2 1 arp 

9oV + V— + -V(p2_)-vx(Vxv) = (3.31) 
2 p op 



In virtue of the identity 



V X (V X v) = V (vV)v (3.32) 



the last equation is equivalent to (1.2). 

Thus, differentiating equations (3.25), (3.26) and eliminating the variables </?, ^, 
77, one obtains the curtailed system (1.1)-(1.3), whereas the system (1.4), (1-12), 
(1.13) follows from the system (3.21)-(3.26) directly (i.e. without differentiating). 
It means that the system (1.4), (1.12), (1.13) is an integrated system, whereas the 
curtailed system (1.1)-(1.3) is not, although formally they have the same order. 

The action of the form (3.20), or close to this form was obtained by some authors 
(Seliger and Whithem, 1967; Salmon, 1988), but the quantities g", a = 1,2,3 are 
always considered as additional dependent variables (but not as indefinite functions 
of I which can be expressed via initial conditions). The action was not considered 
as a functional of fixed indefinite functions g'^{^)- 

The variable rj was introduced, for the action be invariant with respect to the 
transformations of the whole relabeling group (1.16). To understand what the r) 
means from the mathematical viewpoint, let us return to the action (3.9), where 
the internal energy E has the form (3.10). Adding new variables j by means of 
designations (2.16), one obtains instead of (3.14) 

/• ;2 a J 

Ae[P:L^:P]= V—-pE-pk\j''-mpo{i)-—]}dtdyi (3.33) 

where E has the form (3.10). 

Variation of (3.33) with respect to leads to the equation 

^ dE{p,So{^)) dSo ^ .oo4^ 

^ocP^P ^ q; = 1,2,3 (3.34) 

where linear operators are defined by (3.15). Equations (3.34) are linear non- 
uniform equations for the variables p. A solution of (3.34) is a sum of the general 
solution (3.19) of the uniform equations (3.15) and of a particular solution the non- 
uniform equations (3.34). This particular solution depends on the form of the func- 
tion 5*0 and cannot be found in a general form. Adding an extraterm —r]j^dkS with 
7] satisfying (3.24) to (3.13), a reduction of non-uniform equations (3.34) to uniform 
equations (3.15) appears to be possible. Thus, the extravariable r] is responsible for 
the particular solution of (3.34). 

Prom the viewpoint of the action (3.33) a dependence of the internal energy 
E on the entropy means simply a dependence of E on the labels | via a function 
S{i). If such a dependence cannot be expressed through one function (for instance 
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E — E[p, S\{i), S2{i)]) the ideal fluid is described by two entropies Si and S2 and 
by two temperatures Ti = dE/dSi, T2 — dE/dS2- Such a situation may appear 
for a conducting fluid in a strong magnetic fleld, where there are two temperatures 
- longitudinal and transversal. 

Thus flve equations (1.4), (1.12), (1.13) with S", p and v, deflned respectively by 
(1.9), (1.10) and (1.11), constitute the flfth order system for flve dependent variables 
C — {Co,^}, V- Equations (1.1), (1.3), (1.4), (1.12), (1.13) constitute the seventh order 
system for seven variables p, ^, </?, 77, ^S". 



4 Initial and Boundary Conditions 

Boundary conditions describing vessel walls can be taken into account by means of 
a proper choice of the internal energy E[x, p, S) which can include the energy of the 
fluid in an external potential U. 

E^Eo{p,S) + Uit,x), (4.1) 



where U is some given external potential. For instance, let the fluid move inside a 
volume V. Then 

J 0, inside V 

[ 00, outside V 



C/(x) 



Such a choice of the energy E provides that the fluid does not escape the volume V. 

Let us consider the case, when the fluid flow is considered in the space-time 
region fl deflned by inequalities 

Q: t>0, x^>0 (4.2) 

The region Q has two boundaries: I deflned by the relations t — 0, > 0, and 
B deflned by the relations x^ — 0, t > 0. The initial conditions for the system of 
equations (1.1)-(1.4) have the form 

p(0,x)=pi„(x), v°(0,x)=<(x), a = 1,2,3 (4.3) 

5(0,x) = 5i,(x), ea(0,x)=ein(x), « = 1, 2, 3 (4.4) 

at X G X (t = 0, > 0). Here pin, Vin, S'ln, are given functions of argument x. 
The boundary conditions on the boundary B oi fl have the form: 

p(x)|^3=o = Pb(i,y), S(x)ls=,^S^{t,y), {t,y}eB (4.5) 

t;"(x)|,3=o = <(i,y), a = 1,2,3, {t,y} e B (4.6) 

U^)\,s=o^&^y)^ a = 1,2, 3, {t,y}eB (4.7) 

where 

y = {x\x^} (4.8) 
Here pb, Sy,, Vb, a-re given functions of the argument {t,y}. 
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Let us show that indefinite functions g, 5*0, po can be expressed via initial and 
boundary conditions (4.3)-(4.7). The initial conditions for the system (3.21)-(3.26) 
have the form 

ea(0,x)=er„(x), a = 1,2, 3 (4.9) 

P(0,X) = pin(x), ^(0,X) = So[i,n{^)], (4.10) 

(/j(0,x) = (/?i„(x), ?7(0,x) = 77in(x), (4.11) 

(4.9)-(4.11) take place at x e T. The functions </?in(x), 77in(x) as well g°'{^) are to be 
determined from the relations 

= <(x), a = 1,2,3; x e X. (4.12) 

It is clear that five functions g,(fiin,Vin cannot be determined unambiguously from 
three relations (4.12). 

There are at least two different approaches to determination of functions ^in(x) 
and g(^). 

(1) One fixes the functions ^in(x) in some conventional way, sets 

<^in(x) = 0, r/in(x) = 0, xeX (4.13) 

and determines functions g from three relations (4.12). 

(2) Functions g are fixed in some conventional way, and remaining functions are 
determined from relations (4.12) 

The first way. Let the condition (4.9) be given in the form 

ea(0,x)=ein(x)=x", «= 1,2,3, xeX. (4.14) 

In other words, ait — the labels $ coincide with the Eulerian coordinates to within 
a constant factor. The relations (4.12) take the form 

/[Ci,(x)] = v^i^), a = 1,2,3- X e J, (4.15) 

which are resolved in the form 

9"{i) = vm, a = 1,2, 3, e3>0, (4.16) 

Thus, the functions g are expressed through initial conditions (4.3). 

The boundary conditions for the system of equations (3.21)-(3.26) have the form 

ea(a:)L3=o = eb(^,y), a =1,2,3, {t,y} e B (4.17) 

5(a;)L3=o = So[Ut,y)] = ^b(t,y), {t,y} G B, (4.18) 

p(x)|^3=o= Pb(:r) 1^3=0, v(x)|^3=o = Vb(t,y), {ty} ^ B, (4.19) 

(^(x) 1,3=0 = ^(^) 1.3=0 = 0, {t,y} e B, (4.20) 
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Let us set 

C^{t,y)^x", a = 1,2; ^ -ct, {t,y) e B, (4.21) 

where c is a constant. 

Writing relations (1.4) and (3.26) for ,^3 < on the boundary B and using (4.20), 
(4.21), one obtains constraints for the functions g(^) 

/[€b(t,y)]a„e^(t,y) = <(t,y), a -1,2, {t,y}eB (4.22) 

9^[ib{t,y)]do^^{t,y) = -Mt,y), {t,y} e B, (4.23) 

where 

^_ Vb(^,y) I d{f^{t,y)E[p^{t,y),S^it,y)]} r , . ^ 

i^b(t,y) = ^— + , {t,y}eB, (4.24) 

Substituting relations (4.21) into (4.22), (4.23), one obtains three equations for 
determination of functions g(^). Resolving this system of equations with respect to 
g, one obtains 

^"(0=<(-Wc,ei,6), a =1,2; e3<0 

g%^) = c-iXb(-6/c,ei,6), 6 < (4.25) 

Thus, g(^) is determined by (4.16) for ^3 > and by (4.25) for ^3 < 0. In other 
words, the boundary conditions and the initial conditions determine the vector field 
g(^) in different regions of the argument ^. All information about the velocities at the 
initial moment and on the boundary has been transferred into dynamic equations. 
The field g(^) can describe both initial and boundary conditions. 

The second way. Let us choose the functions g in a simple form. Let for instance, 

9\i)=^2, /(0 = 0, 9'{i)^0 (4.26) 

Let us set 

X = V, X = ^2, l^ = ^i (4.27) 
Then the expression (1.11) takes the form 

u(x, A, II, rj, S) = Vx + AV// -rjVS^v (4.28) 

where x, A,/x, are Clebsch potentials (Clebsch, 1857; 1859). Now six equations (1.1), 
(1.3), (3.22)-(3.26), (3.27) [(3.22) for a = 3 is of no importance] for six dependent 
variables p, x, A, /i, rj, S do not contain indefinite functions and have an unambiguous 
form. 

dop + V(pu) = 0, doX + (uV)A = 

doi^ + (uV)At = 0, doS + (uV)S = (4.29) 

dov + (uV)r^ = doX + Xdoii - rid^S + ^u^ + = 
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where u is defined by (4.28). 

The initial conditions for variables p, x, A, /i, r), S are determined by relations 

p(0, x) = pin(0, x), S{0, x) = 5in(0, x), (4.30) 

VXin + AinVpin " 77inVS'i„ = Vin (4.31) 

Three equations (4.30), (4.31) do not determine the initial conditions 

X(0, x) = Xin(x), A(0, x) = Ai,(x), (4.32) 
p(0, x) = pin(x), 77(0, x) = r;i„(x), (4.33) 

unambiguously. 

If the fluid is described in terms of Clebsch potentials, the dynamic equations 
contain neither arbitrary functions, nor information about the initial conditions. 
It should be interpreted in the sense that the description (4.28)-(4.29) in terms of 
the Clebsch potentials is a result of a change of variables in dynamic equations 
(1.1)-(1.3), whereas the description (3.21)-(3.26) is a result of integration of the 
dynamic equations (1.1)-(1.4). In other words, the description (4.28)-(4.29) in terms 
of Clebsch potentials relates to the description (3.21)-(3.26) in the same way, as a 
particular solution of a system of differential equations relates to a general solution 
of the same system. 

Let us note that there are many other ways for determination of indeflnite func- 
tions g(^). For instance, for slightly rotational flows the functions g(|) may be 
chosen as small corrections to the basic irrotational flow described by the potential 
ip. 



5 Description in Lagrangian coordinates 

To show that the system (1.4), (1.9)-(1.13) is indeed the integrated system of hy- 
drodynamic equations, let us rewrite it in Lagrangian coordinates, when flve vari- 
ables X — {x^, x^, x^}, (f, T] are considered as functions of four independent variables 
^ = {^1, ^2, ^3, }■ It is necessary to introduce designations 

Q^ d{x\x^,p d{^^^^ 11^ «'/5 = 1.2,3 (5.1) 

^(^1,6,^3) d{^) dip 

do 

Xa,,^^, «,/3 = l,2,3 (5.2) 



It follows from (2.3), (3.4) that 



Q=(|^r' = ^. (5.3) 



Identities (2.7), (2.8) take the form 

d 



x"'^X^^p = 6^Q, ^^a,7 = 0, a, 13=1, 2, 3. (5.4) 
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Derivative with respect can be recalculated into derivative with respect to as 
follows 



d{^i,^2,^'i) d{x^,x'^,x^) ^ '^^''^9^^ 
or for a derivative with respect to a;" 



^^^g-^X,,^^ a = 1,2, 3 (5.6) 



Applying the rule (5.6) to (1.11), one obtains 

= v^{t. i) = Q'^^o.A^^ + /U) - ^^1^]' « = 1, 2, 3 (5.7) 
Let D/ Dt means derivative with respect to t at constant 4 

^.S„,^^S., (5.8) 

In particular 

^ = aoe„ + ^a^^a = 0, (5.9) 

because D^a/ Dt = by definition. Comparing (5.9) with (1.4) and using (5.7), one 
concludes that 

= ^ = Q-'X^^A^^ + /U) - « = 1> 2, 3 (5.10) 

Now in virtue of (5.6), (5.7), (5.10) the equations (1.12), (1.13) can be rewritten in 
the form 



Dip 1 



o=3 



Dt 2„ , 

a=l 



a,/3 



Liry aE(p, S) 



+ «M = o (5.11) 

9p 



where 

The system of five equations (5.10)-(5.12) is a system for five dependent variables 
X, ?7, considered as functions of four independent variables t, ^. Five indefinite 
functions g,po,5'o of ^ are determined from initial and boundary conditions. In 
particular, if there is no inflow of the fluid and the initial values for x, (/?, 77 are given 
in the form 

x(0,0=^, (5.14) 
(^(0,0 = 0, ?7(0,^) = 0, (5.15) 
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it follows from (5.10), (5.13), that 



S{0, ^) = So{^), p(0, = Po(^), ^(0, ^) = Vin(0 = g(0 (5-16) 

In the Lagrangian description the initial conditions (5.14) for the position of a par- 
ticle labeled by ^ look quite reasonable. From physical point of view a necessity 
of initial conditions for the particle position does not raise doubts. The system of 
hydrodynamic equations (5.10)-(5.12) in the Lagrangian form looks as partly inte- 
grated system. Indeed, the system (1.1)-(1.4) written in the independent Lagrangian 
coordinates ^ for six dependent variables x = x(t, $), v = v(t, $) has the form 

= = -Poi^r'Xa,,Q^^^^[p'E{p, Som, « = 1, 2, 3 

(5.17) 

where X^^/s and p are functions of ^, x"' and x'^'f^ = dx'^/dS,^ which are defined 
by relations (5.1), (5.2), (5.13). The six order system (5.17) contains only two 
indefinite functions po(^) and >S'o(^) describing initial values of density and entropy. 
Initial values of velocity v(0, ^) and position x(0, 4) are given by initial conditions 

x(0,0=^, v(0,^)=Vi„(O (5.18) 

The relation (5.7) is an integral of (5.17). It satisfies the equations (5.17) for any 
functions g in virtue of equations (5.10)-(5.12), although this circumstance is not 
evident directly. 

Note that the curtailed system (1.1)-(1.3) cannot be written in the Lagrangian 
form directly, because it does not contain any reference to Lagrangian coordinates 
^. To introduce ^, it is necessary to append equations (1.4). Then fifth order system 
(1.1)-(1.3) turns to the complete eight order system (1.1)-(1.4). Equations (1.1) and 
(1.3) can be integrated on the basis of (1.4) in the form (1.9), (1.10). The remaining 
equations (1.2), (1.4) constitute the sixth order system which can be written in the 
Lagrangian form (5.17). 



6 Incompressible fluid 

In the special case of the incompressible fiuid, it should set p = po =const in the 
action (3.20) and introduce new variable 

v = j/po, Po = const (6.1) 

It is easy to verify that rj = rj{^), S = S'o(C), and the last term of (3.20) can be 
incorporated in the term j^g°'{i)dk^a- Thus, the action for the incompressible fluid 
looks as follows 

2 

^e[v, ^, 97] = Po /{ y - vV</p - 9''{^)doia - 9''U)vVU}dtdx, (6.2) 
where g°'{^) are arbitrary flxed functions of ^. 
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Variation with respect to v, ^, 93 gives 

5v: v = V(^ + 5"U)Vea (6.3) 

Po = {9'''^ - g^nma + vV{„) =0, a = 1, 2, 3 (6.4) 

Po'^ = Vv = (6.5) 

In the general case the condition (3.27) is satisfied. Substituting (6.3) into (6.4) 
and (6.5), one obtains 

SoCa + [^V + /(^) VC/?] VCa = 0, a = 1, 2, 3 (6.6) 

VV + (f'^m^p^ia + ^?"(^)V2a = (6.7) 

The dynamic equation for 93 does not contain temporal derivative. 
Conventional hydrodynamic equations for the incompressible fluid 

Vv = 0, 9oV + (vV)v = (6.8) 

Po 

are obtained from relations (6.3)-(6.5). Differentiating (6.3) with respect to one 
obtains 

aov = V[ao^ + g'\i)doia\ - fi^^^oOV^a (6.9) 
where fi"'^ is defined by (1.21). It follows from (6.3), (3.28), (3.29) 

vx(V X v) = 0"^(4)Ve/3(vV)ea. (6.10) 

In virtue of (6.4) the last term in rhs of (6.9) coincides with rhs of (6.10). Then 
using the identity (3.32), one obtains 

+ (vV)v = V[ao(p + ^"(O^oCa + ^v^] (6.11) 



The equation (6.11) coincides with the second equation (6.8), provided one uses 
designation 

V-^o<p-^"(^)9oea (6.12) 



P 1, 2 



Po 2 

Here the pressure p is determined after solution of the system of hydrodynamic 
equations (6.3)-(6.5), or (6.8). 

If the condition (3.27) is fulfilled, the equation (6.4) may be written in the form 
(6.6). Then eliminating 9oCa by means of (6.6) and using the corollary of (6.3) 

5"(^)(vV)ea = v2-(vV)v3, (6.13) 

the relation (6.12) for the pressure can be written in the form 

-v2 - [^09? + (vV)<^] (6.14) 



P 1_ 2 



Po 
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In the case of a irrotational flow, wfien v = the inequahty (3.27) turns 
to equahty and formally the derivation of the expression (6.13) is not founded. 
Nevertheless the relation (6.14) remains vahd in this case, because it turns to the 
integral 

— + + doif = (6.15) 
Po 2 

The relation (6.14) is vahd for any flow of incompressible fluid, but it is not an 
integral. It is a definition. Prom point of view of the description in terms of hy- 
drodynamic potentials (/?, |, the relation (6.14) is a definition of the pressure p in 
terms of dependent variables (/?, | and the relation (6.3). From the viewpoint of the 
curtailed system (6.8) the relation (6.14) is a definition of the function 99 in terms 
of variables v, p. Nevertheless, this definition is useful for a description of a slightly 
rotational fiow, when the velocity v can be represented in the form 

v = V<^ + 5v, |5v|<|V<^|, 5v = 5"(^)Vea, (6.16) 

where 6\ describes a small rotational component of the velocity. 
Substituting (6.12) into (6.14), one obtains 

^^-do^-\{V^f + \{5^f (6.17) 
Po 2 2 

This relation shows that a contribution of the small rotational component 5v into 
the pressure is of the second order (5v)^ and this contribution always increases the 
pressure. 



7 Slightly Rotational Flow of Incompressible Fluid 

It seems to be reasonable to consider a slightly rotational fiow as a small correction to 
an eatablished irrotational fiow which can be effectively calculated for a fiow around 
different bodies. Let m be a set of dependent variables, and dependent variables uq 
describe an irrotational fiow. Let us represent a slightly rotational fiow in the form 

u^uq^sui^ 0{e^) , £ < 1 (7.1) 

where e is a small formal parameter which is set to be equal to 1 after calculation. 
It means that the quantities Ui are considered as small with respect to Mq, Ui/uq = 
0{e) <S 1. Substituting (7.1) into dynamic equations and neglecting higher order 
terms, one obtains linear equations for Ui. Use of dependent variables v, p, satisfying 
equations (6.8), leads to the following first order approximation equations. 

Vvi = 0, 9oVi + (V(/?oV)vi + (viV)V(^o = -Vpi/po, Vq = V^o (7.2) 

Dynamic equation for the rotational component wi = V x Vi of the velocity follows 
from the second equation (7.2) 

doui + (V(/?oV)a;i - (wiV)V<^o = (7.3) 
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Both equations (7.2) and (7.3) are difficult for a solution. Use of dynamic equations 
(6.3), (6.6) and (6.7) for hydrodynamic potentials appears to be more effective, 
because for a fixed velocity v = V<^o equations (1.4) for ^ are equivalent to the 
system of ordinary equations (1.6) which can be solved simply enough. Indeed, let 
the velocity v = V<^o of the basic established irrotational fiow be known, and the 
orthogonal coordinate system C,r),C coupled with the flow be determined by the 
relations 

e = «-Vo(x), Ve-V77 = 0, Vr7-VC = 0, • VC = (7.4) 

where u^c is a constant velocity of the flow at inflnity. Then the general solution of 
the equations (1.4) is written in the form 

fa{X -u^t,riX): a = 1,2, 3 (7.5) 
X = X{C, V,0-J Vo-'ie, V, Ode, = (VO' = (VW^oo)' (7.6) 



where f^, a = 1,2,3 are arbitrary functions of three arguments X — Ur^t, 7], 
Vo = Vo{^, 1], () is a dimensionless velocity of the basic fiow considered as a function 
of coordinates {C,vX)- This fact can be tested by means of a direct substitution 
of (7.6) into equation (1.4) with v = Vcpo. X can be interpreted as a "distorted" 
Cartesian coordinate in the direction along the basic flow at inflnity. The distortion 
of the coordinate X is chosen in such a way that the fiow were uniform. 
Let us use the expansion (7.1) for dependent variables u — {ip, ^} 

V = Vo + £Vi + Oie^) = V ((/^o + sifi) + eg^i^)V^a + 0(£') (7.7) 

Vo = Vv^o, vi = Vv?i + (?°(^) (7.8) 

Substituting (7.7) into (6.7), (6.6) and equating coefficients before equal powers 
of e to zero, one obtains 

£0 : vVo = (7.9) 

e: ^oCa + Vv^oVa = 0, a = 1,2,3 (7.10) 

e: V2a + V[(7"(OVa] = 0, (7.11) 

Let us note that the labels | are calculated only in zeroth approximation. Instead of 
the expansion of ^ one uses the supposition that the rotational components g"'{^) of 
the velocity are small as compared with the irrotational component V</?o- General 
solution of equations (7.10) have the form (7.5). Then (7.11) takes the form 

VVi + V (A^VB^) = 0, ^"(^)VC„ = A^VB^ (7.12) 

where Aa, Ba, a = 1,2,3 are arbitrary functions of arguments X — Uaot, i], (. Let 
S be the solid boundary of the volume with the fiuid. Then equation (7.12) should 
be solved with the boundary condition 

nV(fii I s = - Ac,n VSc, I s (7.13) 
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Note that 991, ^ are hydrodynamic potentials which are not determined inam- 
biguously by the fluid flow. But any set of hydrodynamic potentials determines 
inambiguously a fluid flow. It is possible one to express the functions A = {^a}, 
B = {Ba}, « = 1,2, 3 via functions g(^) — ct = 1,2,3, choosing at t = 

^ = X. Then the relations 

B(X,77,C)=x, A(X,77,C) = g(x) (7.14) 

determine the form of functions A, B and the fact that A is a function of B: 
A=g(B). On the other hand the functions g(^) are indeflnite functions whose 
form is determined by the initial conditions and by a choice of the labeling ^. The 
relabeling transformation (1.16) changes the form of the functions g(C)- 

Besides only the rotational component of the vector AaVBa is essential, because 
its irrotational component is compensated by the contribution of (pi into v. Indeed, 
Let AaV BaCan be represented in the form 

A^VBa = V$ + V X C, VC - (7.15) 

where C is some solenoidal vector describing the rotational component of AaVB^- 
Substitution of (7.15) into (7.12), (7.13) leads to the relations 

VVi = 0, = + nV(^i|s = - nV X C|s (7-16) 

In other words, the potential Lpi compensates the irrotational component of 
AaVBa- As a result the total potential 

tildevarphii satisfles the same equation (7.9) as (po does. The boundary condition 
for (fi appears to be qenerated by the rotational component of Vi 

v = V{(po + e(fi)+V xC + 0{e^), V\ipo + eifi) = (7.17) 

nW{(pQ + e(fi)\.^^ - enxC\j. (7.18) 

Solution of the equation (7.12) with the boundary condition (7.13) is written in the 
form: 

<^i(x) = ^1 G(x,x')V'K(x')VX(x')]rfx'- ^ ^G'(x,x')A„(x')V'5„(x')rfS' 

(7.19) 

where V means the gradient with respect to coordinates x', and G(x, x') is the 
Green function satisfying the following conditions 



7/2r^. n_ M aG'(x,xO 



V''G(x,x') = -47r(5(x-x'), 



dn' 



= (7.20) 



where d/dn means component of the gradient in the direction of the normal, and 
functions Aa, Ba depend only on arguments X—Uoot, r), C which are known functions 
of X. The surface integral is taken over the surface surrounding the volume V 
with the fluid. This surface includes the surface S of the body and the inflnite 
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surface around V. The surface integral in (7.19) contains functions Aa{'x.'), Baipd) 
determined in the whole volume V . It permits to transform the surface integral to 
the volume one 

y G'(x,x')A«(x')V'5«(x')(iS' = V'[G'(x,xOA«(x')V'5„(xO](ix' (7.21) 

V 

Then the relation (7.19) transforms to 

(^i(x) = / V'G'(x, x') • A«(x') V'S„(x')dx' (7.22) 

V 

Thus, the shghtly rotational flow is described by the relations 

V = V<^o + £ (V(^i + A^VB^) + O {e^) (7.23) 

^ = -9o (^0 + e^i)-]- ( Vv^o)' - £V<^o V<^i + O ie'') (7.24) 
Po ^ ^ 

where A^, a = 1,2,3 are arbitrary functions of arguments X — u^t^ 1], ( and 

ipi is described by the relation (7.22). Expressions for p is a result of application of 

(6.17) to the case of the slightly rotational flow. 

Relations (7.23), (7.24) satisfy the equations (7.2), (7.3), but a direct deriva- 
tion of (7.23), (7.24) from (7.2), (7.3) seems to be rather difficult because of terms 
(a;iV)V(/9o! (viV)V(/9o- Such terms connect three similar equations of the type (1.4) 
into united system of three equations and prevent to use equivalence of one equation 
of the type (1.4) and the system (1.6) of ordinary equations. 

The quantities ^q(^), Ba{^) and other functions of only ^, for instance Q"^^(^), 
defined by (1.21), depend on only arguments X — Uoc,t, 1], ( and can be regarded as 
quantities frozen in the basic flow, because they depend on t and X only through 
X — Uoot- Gradients of frozen quantities, for instance AaVBa, are not frozen in 
the basic flow, in general, because they depend on t and X not only via argument 
X — Uoot. They depend on X also through arguments • dX/d^ — V^/ |V^|^, Vr^, 
VC which depend, in general on X, but do not depend on t. It is a reason why a 
description in terms of "frozen quantities" appears to be more effective, than that 
in terms of the "frozen quantities gradients". Appearance of extraterms in (7.2), 
(7.3) as compared with (1.4) is connected with the use of the gradients of frozen 
quantities. 

8 Wave Function and Kinematic Spin of a Flow 

There exists a special complex form of hydrodynamic potentials. In this form the 
dynamic equation for the irrotational component of the flow is very close to a linear 
equation. 

Idea of the transformation of hydrodynamic potentials is very simple. Let there 
be an irrotational flow v = 6V0 of a compressible fluid whose internal energy de- 
pends only on p. Here 6 is a constant of dimensionality [6] = [L^T~^] which is 
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introduced to make the potential (f) to be dimensionless. Then dynamic equations 
(1.1), (1.2) can be integrated and written in the form 



7 2 a 

dop + bV{pV<l>)^0, bdo<j> + -{V<py^- — [pE{p)] (8.1) 

Equation (1.2) is obtained as a gradient of the second equation (8.1). The term 
6^(V0)^/2 is the principal nonlinear term known as convective nonlinearity term 
(vV)v. Introducing the complex variable ^, this term can be removed. Let 

* = ^e^^ p = v = -^(W*-V**) (8.2) 

The bar over symbol means the complex conjugation. Then equations (8.1) are 
described in the form of one complex equation for the complex dependent variable 

zbdo^ + |v^vl/ = {^[pE(p)] + -^v2(Vp)}^, P^^^ (8.3) 

Equations (8.1) are imaginary and real components of (8.3). The second term in 
Ihs of (8.3) is linear. It corresponds to the convective term (vV)v. Now this term 
contains the second order spatial derivative. It is a price which is paid for the 
linearity. As a result of the transformation (8.2) the nonlinearity is transmitted from 
the kinematic term to the dynamic one. The kinematic nonlinearity is stronger, 
than the dynamic nonlinearity, because the kinematic non-linearity is connected 
directly with velocities, whereas the dynamic nonlinearity is coupled directly only 
with forces and with accelerations. Thus, a displacement of the nonlinearity from 
kinematic terms to dynamic ones weakens the nonlinearity (Rylov, 1989). 

The complex variable is used in quantum mechanics, where it is known as 
the wave function. The fact that the Schrodinger equation for the wave function 
describes a irrotational flow of some ideal fluid is well known (Madelung, 1926). For 
such a fluid the internal energy depends on p and Vp. This internal energy has such 
a form {E = 6^(Vp)^/2p^, b = —h/m) that rhs of (8.3) vanishes, and the equation 
(8.3) becomes linear. 

In the case of incompressible fluid rhs of (8.3) becomes indeflnite due to additional 
constraint. The irrotational flow of the incompressible fluid is described by the 
equations 

a§ + f V^* = iL*. (8.4) 

at 2 po 

f ^ = Po = const (8.5) 

where the pressure p is considered as some function of time and position. The 
situation is the same as in the case of equations (6.8), where the value of the pressure 
p is determined only after the flow is calculated. 

The equation (8.4) looks as a linear equation. But in reality it is nonlinear 
because of the nonlinear constraint (8.5). In the case of a rotational flow some 
additional nonlinear terms appear. However if the flow is slightly rotational, there 
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is a hope that these additional nonhnear terms, as well the nonlinear constraint 
(8.5) could be considered as corrections to the irrotational flow described by linear 

dynamic equation. 

Using idea of linearization of the irrotational flow, let us introduce n-component 
complex function ^' = {^'a}, « = 1, 2, ... n, defining it by the relations 

= VPe^'^/VU), ^a = Vpe-'^/'Wo,{^), t*=X:*a*a (8.6) 

a=l 

where the bar over a symbol means the complex conjugate, Wa{^), a = 1, 2, ... n are 
complex functions of only variables ^, satisfying the relations 

-|E(*a^-^^.)=/(^), /3 = 1,2,3, tw^w^^l (8.7) 

n is such a natural number that equations (8.7) admit a solution. Here b is some 
constant of dimensionality [b] = [L^T^^] which is introduced to make the quantities 
Wa and (p = (fi/b dimensionless. The number n and the form of functions w depend on 
the form of functions g". But if Wa, a = 1, 2, ... n is a solution of (8.7), the variables 
\1/q, are some functions of p, ip, ^, but not of their derivatives, although equations 
(8.7), determining the transformation (8.6), contain derivatives of functions Wa with 
respect to |. 

It is easy to verify that p and pv with v defined by (1.11) have the form 

ib - - 

p = j = pv = --(W*-V*^)-77V,5** (8.8) 

Then the variational problem with the action (3.20) appears to be equivalent to the 
variational problem with the action functional 

/ib - - - 

{-{Wo^ - do^ ■ ^) + ridoS^^i/ 

1 ib - 

[-(*V* - • ^) + r]VS ■ - E[x, (8.9) 



2** '2 

Note that the function considered as a function of independent variables {t, x} 
is very indefinite in the sense that the same fiuid fiow may be described by different 
\l'-functions. There are two reasons for such an indefiniteness. First, the functions 
Wa{^) are not determined uniquely by differential equations (8.7). Second, their 
arguments ^ as functions of x are determined only to within the transformation 
(1.16). Description of a ffuid in terms of the function \& is more indefinite, than the 
description in terms of the hydrodynamic potentials ^. Information about initial 
and boundary conditions containing in the functions g(^) is lost at the description 
in terms of the ^'-function. The two-component ^-function can be obtained directly 
from the Clebsch variables by means of a proper change of variables (Rylov, 1989). 
However, one cannot be sure that any fiow can be described by two component wave 
function. According to (8.7) the classification of wave functions over the minimal 
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number of components is connected with the form of functions g(^) and, hence with 
the integration of hydrodynamic equations. 

Let the function \I/ have n components. Regrouping components of the function 
^ in the action (8.9), one obtains the action in the form 

+ T E QaP,^QaP,^P + ^(Vp)' - pE}d\ p = (8.10) 

4 a,f3=l 



where 
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^ = {^,A}, Ao = 7)8^3, A = riVS, (8.11) 
^ ^" a,/3= l,2,...n 7 = 1,2,3 (8.12) 



Corresponding dynamic equations have the form 

5A 1 6^ n _ 

^ u=l ^ u=l C'p Sp 

-V(^Vp)*„ = 0, a = l,2,...n (8.13) 

^ = -9,(/5) = 0, (8.15) 

where j = {p, j} = {j''}, /c = 0, 1, 2, 3 is defined by (8.8). 

In the case of the irrotational flow, when g°'{$) — d^{^)/d^a equations (8.7) 
have a solution for n = 1, and the function ^ may have one component. Then all 
Qai3,'y = 0, as it follows from (8.12). 

The number n of the ^-function components in the actions (8.9) and (8.10) 
is arbitrary. A formal variation of the action with respect to and ct = 

1, 2, . . .n leads to 2n real dynamic equaitons, but not all of them are independent. 
There are such combinations of variations 5^q,, S^a, ct — 1,2,. . .n which do not 
change expressions (8.8) and ^(?o^ — ■ \E'. Such combinations of variations 5^q,, 
S^a, OL = 1, 2, ... n do not change the action (8.9), and corresponding combinations 
of dynamic equations bAjb'^Q, = 0, bAjb'^Q, = are identities that associates with 
a dependence between dynamic equations. 

Thus, increasing the number n, one increases the number of dynamic equations, 
but the number of independent dynamic equations remains the same. In such a 
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situation it is important to determine the minimal number rim oi the ^-function 
components, sufficient for a solution of equations (8.7) with the given vector field 

g^{^) in the space of the labels |. 

Note that under the relabeling transformations (1.16), the quantity g(^) trans- 
forms as a vector 

/(O-^^® = §^^"U), /3 = 1,2,3 (8.16) 

It is necessary for the quantities (8.8) and the action (8.9) to be invariant with 

respect to the transformation (1.16) 

Let ^ be a set of all vector fields g^{i) in V^, and Qn be a set of such vector fields 
g^{^) in which can be represented in the form 

n 

9^^) = j:vl{mk{^)/dC(S: Vi^l (8.17) 

k=l 

where n is a fixed natural number, and the functions rjk, (k, k = 1,2, . . .n are scalars 
in V^. Under the relabeling transformation (1.16) the functions (8.17) transform as 
follows 

Vk{^)^'nk{l)^Vk{^), a(0 ^ Cfe(^) = CfeU), k^l,2,...n 

A^) - m) = §^5"U) = ^tvli^)^ = tm^ (8.18) 

In other words, a vector field g^{$,) of the form (8.17) transforms into the vector 
field g'^(^) of the same form (8.17), and the set Qn is invariant with respect to the 
group (1.16) of the relabeling transformations. 
It is easy to see that 



^Qn, 6^0 = 0, n=l,2,... (8.19) 



because the nth term of the sum (8.17) can be combined with the first one, if C„ is 
a function of rjn- Let 

Sn = Gn\Gn-l, 71=1,2,... (8.20) 

Then _ 

S — Tim 

g= [j Ss, Si = ^, / = n^ + l,n^ + 2,... (8.21) 

s=l 

where rim is the number of nonempty invariant subsets of the set Q. Each subset 
Sk contains only such vector fields g^{^) which associate with the A;-component 
function ^ = {^a}, a = 1, 2, ... /c, having the components 

*i = |(l-E^a)p| exp[i(0 + Ci)], 
*a = ?7av^exp[i(0 + Ca + Ci)], « = 2, 3, . . . A; (8.22) 



30 



In particular, the set Si associates with an irrotational flow, described by a one- 
component ^-function determined by one scalar Ci; and the set S2 associates with 
a rotational flow described by a two-component ^-function, determined by three 
scalar functions (i, 772, C2 (Clebsch variables). 

Thus, types of the perfect fluid flows can be labeled by invariants of the relabeling 
group (1.16). This labeling is connected with the minimal number of the wave 
function components. In the quantum mechanics the minimal number of the 
wave function components is connected with the spin s of a particle described by 
this wave function. The relation has the form 

s^ium- l)/2 (8.23) 

We shall refer to the quantity s deflned by (8.23) as the kinematic spin (k-spin) 
of the fluid flow. According to this deflnition the kinematic spin of an irrotational 
flow is equal to 0. The kinematic spin of the rotational flow s > 1/2. The k-spin 
s = 1/2 is possible. Is a flow of higher k-spin possible? This interesting question is 
yet open. It seems to be connected with the complicated problem of knottedness of 
vortex lines (Moffat, 1969; Bretherton, 1970). 

As far as the fleld g(^) is determined by the velocity held Vin(x), at the initial 
moment, when ^ = x and g(x) = Vin(x), the k-spin is determined by the initial 
velocity Vin(x). In particular, if the initial velocity can be represented in the Clebsch 
form (4.31) with rji^ = 0, the k-spin of such a flow is equal to 1/2. 

Let us consider a description of the incompressible fluid in terms of the wave 
function. For the incompressible fluid the action (8.9) has the form 

A[^, ■^,p]^poll {'-^{^0-^ - do^ ■ *) - ^v2+P(l - ^■^)}dtdx, (8.24) 
where p = po =const is the fluid density, 

v = --(W*- •*) (8.25) 

is the fluid velocity, and P is the Lagrange multiplier which introduces the constraint 

p/p, = = 1. (8.26) 

The dynamic equations arc obtained as a result of a variation with respect to P, ^, 
^. Variation of (8.24) with respect to P gives (8.26), and 

p-i^ = i69o* + i^[vV* + V(v*)] - P* = (8.27) 

Let us convolute ^ with (8.27) and take imaginary and real parts of the obtained 
relation. In virtue of (8.26) one obtains respectively 

Vv = --V(*V* - • *) = (8.28) 

P = i6(*9o* - • *) - (8.29) 
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Then the dynamic equation (8.27) takes the form 

ib{do<i/ + vV*) - = (8.30) 

where v is defined by (8.25). Let us use the expressions (8.6), (8.7) of the wave 
function via functions g(^). Using p — 1, one derives 

P = -do^-g''U)do^a-^' (8.31) 

Comparing (8.31) with (6.12), one derives the relation between P and the pressure 
P 

^ = P + vV2 (8.32) 
Po 

9 Perturbation Theory for the Flow of /c-spin 1/2 

Let us consider the flow of k-spin 1/2, when the wave function has the form 

* = (^,7) (9.1) 

If the wave function is one-component, k-spin is equal to and the flow is irrotational. 
If 7 ^ /9, the flow is shghtly rotational. The expanded form of (8.30) 

ib{do(3 + vV/3) - P/3 = 0, ^6(^07 + VV7) - P7 = (9.2) 

v = -^(;gV/3 + 7V7- V;9-/3- V7-7) (9.3) 
Let the expansion of dynamic variables /3, 7, v, P have the form 

/? = /?o + e/9i + e'/?2 + O(e^), 7 = ^Ti + + ^(e^), 

P = Po + ePi + e2p2 + 0(e^), v = Vq + evi + eVa + O(e^), e<l (9.4) 

where e is a formal small parameter. The expansion is used to within e^, because 
it appears that the expansion of variables /3, v, P contains only even powers of e, 
whereas the expansion of 7 contains only odd powers of e. 

Let us substitute expansions (9.4) into dynamic equations (8.26), (8.28), (8.30) 
and equate to zero coefficients before different powers of e. In the zeroth order 
approximation one derives 

;3o/3o = l, /9o = e*^°, Vq = 6V0o (9.5) 

where 0o is a real variable. Then the equations (8.28) and (8.30) give respectively 

VVo = 0, Po = -6^000 -&'(V0o)', (9.6) 

Thus in the zeroth approximation we have a problem of an irrotational incompress- 
ible fiow. This problem is supposed to be solved. 
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In the first order approximation we have 

/9o/9i + /9i/9o = /?i = i0i/?o (9-7) 
where 0i is a real quantity. Substituting (9.7) into (8.28), (9.2), one derives 

V'(/)i = 0, Pi = -69001 - 6'V0oV0i, (9.8) 

As far as the dynamic equation for 0i has the same form as for (po, one can replace 
00 + v^0i, by one variable 0o- It means that formally one may set without a loss of 
generality 

01 = 0, Pi = 0, /5i = (9.9) 
The second equation (9.2) reduces to the form 

i(9o7i + 6V0oV7i) + [9o0o + fe(V0o)']7i = (9.10) 

It is easy to see that (9.10) reduces to the equation of the type of (7.10). The general 
solution of (9.10) has the form 

7i ^ Ri{X - Uoot, r}, C) exp[i0o + i-^iiX - Uoot, r}, ()] (9-11) 

where Ri, 'j?i are arbitrary functions of arguments X — u^ot, rj, ( defined by relations 
(7.4)-(7.6) with cpo = 60o- 

Taking into account (9.9), one has in the second order approximation 

^o/?2 + /32/3o + 7i7i = (9.12) 

This equation has the general solution 

/?2 = (-^7i7i + #2)/3o (9.13) 
where 02 is a real variable. Equation for 72 has the form 

i(ao72 + 6V0OV72) + [^000 + 6(V0o)^]72 = (9.14) 

which coincides with the equation (9.10) for 71. It means the £72 may be included 
in 7i and without a loss of generality one may set 72 ~ 0. 
The equation (8.25) and (8.28) are written in the form 

V2 = 6V02 + bRjV^i, V^02 = -V(P?Vi?i) (9.15) 

Then the first equation (9.2) gives 

P2 = -6(^002 + VoV02 + V2V0o) (9.16) 

From (9.5), (9.15) one derives 
v^V{<po + e^ip2)+e%RlV^i + 0{e^), ¥^0 = &0o, </'2 = &02 (9.17) 
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Using (8.32), (9.6), (9.8), (9.9), (9.16), (9.17), one derives after calculations 



= -do{ipo + eV2) - T^i^^of - e^VipoVip2 + O(e^), (9.18) 
Po ^ 

Results (9.17), (9.18) agree with (7.23), (7.24), provided one sets = e. The expres- 
sion (7.23) contains the combination AaVBa, a = 1,2, 3, whereas the expression 
(9.17) contains only two arbitrary functions i^i. The fact is that that (7.23) 
describes arbitrary slightly rotational flow, whereas the relation (9.17) describes 
slightly rotational flow of k-spin 1/2. For any flow of k-spin s < 1/2, the six func- 
tions Aa, Ba, a = 1,2,3 reduce to two functions (for instance, Ai, Bi). But the 
question about maximal value of k-spin is open now. On one hand, one can prove 
(see, for instance, appendix to the paper by Eckart, 1960) that any vector field g(x) 
in the three-dimensional space can be represented in the Clebsch form 

g = VCi+r/2VC2 (9.19) 

On the other hand there are examples of vector fields which cannot be represented 
in the form (9.19) inside the whole space (Moffat, 1969), if Ci, ''72, C2 are considered 
as single-valued functions of coordinates. 



10 Two-dimensional irrotational flow 

Let us consider a fiow around a circular cylinder of a radius a. The cylinder axis 
is supposed to coincide with the axis of z. The fiow is directed along the axis of 
X. The basic irrotational fiow is described by the potential tpo, and by the stream 
function which are connected with the velocity components Vq = {vxiVy^Vz} by 
means of relations 

Vx = dx^pO = dyljj, Vy = dyipO = "^X^', = (id) 

We are going to use variables ip, ip as coordinates, and introduce another variables 
^,f],C defined by (7.4). In the given case they are determined by relations 

C = "w^Vo = bu^(l)o, ri = bu^xjj, (^z (10.2) 

where is the velocity of the fluid at infinity. In absence of a circulation one 
derives for the circular cylinder (see, for instance Lamb, 1932) 

^=(1-1 — ^)rcos'i?, ^ = (1 ^)rsini?, (10.3) 

r"^ — x^ + y"^, •& — arctan — 

X 

dx^ = dyr), dyi^-dxT) V^e^O, V^r^^O, (10.4) 
In virtue of (10.4) 

Ve • Vr^ = 0, (VO' = = (Vvf = = Do (10.5) 
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where Vo = | V^| /mqo is a dimensionless velocity. For the case of the circular cylinder 
the determinant 

T^o' = i5o = l + ^-^cos(2^) (10.6) 

tends to 1, if r — > oo. 

Let us represent the variable X in the form 

X = e + A, A = A(e, v) = J [Do\e, v) - m' (lo-r) 



Then according to (10.6) 

1.. ^ aV^^ - 2 cos(2?9) 

At ^ ^ ±oo 

A(e,77) = A±(77) + 0(r'), e^ioo (10.9) 

Let us try to find a slightly rotational stationary flow. If the flow is stationary, 
the functions A^, in (7.23) do not depend on argument X — u^ot- are linear with 
respect to X — u^yot and do not depend on r}, (. Let the flow do not depend on ( — z. 
Then 

Ai = f{v), /(r/)^0, 77 ^oo ^2 = ^3 = (10.10) 

Bi=X-Uoot, ^2 = ^3 = (10.11) 

where X is determined by (7.6). Then according to (7.23), (10.2), one derives 

V = u^[l + efiv)]VC + £/(r?) VA(e, v) + eV^^ (10.12) 
The potential is determined by the relation (7.22) which takes the form 

^litri) = -^l V'G(x,xO/(r/OV'A(r,r7')cix' (10.13) 

where ^,77 are functions (10.3) of x = {x,y,z}, and ^',77' are functions (10.3) of 
x' — {x', y', z'}. Introducing designation 

and using (10.4), (10.5), the relation (10.13) reduces to 

M^, ^) = / ^'^(e, V, e', 770/(^0 V'A(e', 77') W (10.15) 
where the Green function G{^, rj; rj') has the form 



G(e, m e', V') = -2 In V(e - e)2 + - 2 in V(C - i'f + (r? + V'f (10.16) 
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Such a reduction of the Green function is possible only in the case, when A (^,77) 
does not depend on ( — z. Here ^, rj are functions of {x,y) defined by (10.3), and 
rj' are the same functions of {x', y'). Calculation of Dq and of dA/d^ in terms of 
variables C^vX gives 



ai + \Jl + a\- 2q!2 + y 2(0;^ + cti^l + - 2q!2 - Q!2) 



;i0.17) 



CKi + \j'2.{a\ + ai^l + ctf - 2a2 — 0.2) 3 

9e ^ 



(10.18) 



Using the relations (7.24), (10.2), (10.5), (10.12), (10.15) and the second relation 
(7.20) written in the form 



(e,^;r,^')U = o 



(10.20) 



one obtains for the pressure p 

p 1 



1 e d f.dGdA' dGdA'.^. , 



(10.21) 
(10.22) 



Now let us calculate the force F acting on the cylinder. Components of the 
external normal n are described by relations 



n = {na:,ny} = {cos'j?, sin'j?} = {^/2a, \Jl- ^^/^o?} 
Components of the force F have the form 



(10.23) 



Po{^)nxad^ 



Fy = - j Po{^)nyadi9 = j p{Odi 



Vl-e/4a2' 

(10.24) 

where the contour integral means the integral in the complex plane of ,^ + iry, taken 
around the cut 77 = between the points a and —a. p(^) = po{9) is the pressure on 
the surface of the cyhnder. 

Using expression (10.21) for the pressure p, one derives 



TT 

F^^-p^ul I cos^ad^?|L>o(e,0) 



]+0(s% (10.25) 
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TV 



} + 0(s% (10.26) 



where 



-oo — oo 
,2 



^ = 20 0081? and 1^0(^)0) = 4sin i). The first term in (10.25), (10.26) vanishes and 



An 



—oo — oo 



dlydA' dIydA'\ . , , 2n 



where 



TV 

— TT 

TT 



5G(C,0;r,V) 



= -2 In 



a cos i^rfi? 



a sin -dd-d 



{2acos'&-^'f + r}'^ 



;i0.28) 
;i0.29) 

;io.3o) 

[10.31) 
;i0.32) 



It is easy to verify that Iy{^', r)') — 0. Then Fy — 0. 

^xiCv) is a single-valued function of arguments and rj^ (see details in Ap- 
pendix). As it follows from (10.6), (10.17), (10.18) dA/d^ is a single- valued function 
of arguments ^'^ and 77^, and A has the form 



where i? is a single- valued function of arguments ^'^ and r]^. Then 



(10.33) 



"'"''^ (C.{e.n% '¥^ = (C,((W) (10.34) 



dr] dr] 



where Ci and C2 are single- valued functions of arguments and r]"^. It follows from 
(10.28) and (10.34) that also vanish in the first order approximation F^ = 0{e^). 
In other words, the D'Alamberian paradox takes place also in the case of the slightly 
rotational stationary flow. 



11 Concluding remarks 

Basing on the Jacobian technique and on the invariance with respect to the relabel- 
ing group, one succeeded to integrate dynamic equtions for the ideal fluid. Three 
indefinite functions g(^) of three arguments arise as a result of this integration. 
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These functions can be expressed vias initial and boundary conditions, and all es- 
sential information on the fluid flow appears to be concentrated in the integrated 
system of hydrodynamic equations. 

Although the integrated system and the curtailed Euler system contain the same 
number of equations, they differ in the relation that the integrtated system takes 
into account information contained in the Lin constraints, whereas the Euler system 
ignores it. This difference is of no importance for irrotational flows, when the fluid 
properties described by the Kelvin's theorem (and Lin constraints) are inessential, 
because they are fulfilled automatically. But this difference is essential in the case 
of irrotational flows. 

One should expect that in the case of a compressible fluid the hydrodynamic 
potentials and the wave function will describe effectively an interaction between the 
acoustic waves and the slight fluid vorticity. 

Integration of dynamic equations generates two ways of descriptions: a descrip- 
tion in terms of hydrodynamic potential (DTHP) and a description in terms of wave 
function (DTWF). Both DTHP and DTWF take into account Lin constraints and 
appear to be effective in the case of shghtly rotational flows. There is a hope they 
will appear to be effective also in the case of strongly rotational flows. Maybe, an 
application of DTHP and DTWF will allow to formulate the turbulence problem 
in a proper way. At any rate DTHP and DTWF seem to be more effective for ro- 
tational flows, than the conventional description in terms of velocity, because they 
take into account information contained in Lin constraints (and Kelvin's theorem). 
DTHP can take into account even initial and boundary conditions. 

The integration of hydrodynamic equations and appearance of the field g ac- 
tivates the relabeling group. Invariant subsets of this group can be used for a 
classification of the fluid flows. The field g appears to be a tool for introducing 
such attributes of the ideal fiuid as the wave function and the spin. These concepts 
are new for conventional hydrodynamics, although they are well known in quantum 
mechanics. In some cases it may be useful in the sense that quantum mechanical 
methods can be used in hydrodynamics and vice versa. 

The author is indebted to Prof. V.A.Gorodtsov for fruitful discussions which 
were very useful for writing this paper. 
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12 Appendix 

Let us calculate the expression (10.30) 
UC,V) = J i^o(2acos^,0) 



acosM'd {A.l] 



C'=2acosi? 



where 



G{2a cos i?, 0; r/) = -2 In [(2a cos^ - + rf 

dG _ 1 dG _ 4(2acos^9-0 
5^' 2a sin 1? Si? {2acosd — + 

Let us introduce designations 

2a 2a 



Then 



4(e,^) = 8 



4 y"(l - cos^ 1?) COS 



Do (2a COS ^^,0) = 4sin2't? 
cos?? — fli 



{cos-d — iiiY + 
( 1 



1 - cos ■»?) COS M'^ = 



+ 



1 



cos ■»? — /^i + i//2 COS — III — ijJt,2 

— n 

This integral reduces to the form 



(A2) 
(.1.3) 

(A.4) 

(A5) 



(^.6) 



J y cos v — 11 cos v — 11 

— TT 



where 



A = (//^ + /i^ + ///i - 1) = 3//? - A«2 - 1 
Integration of (A. 7) leads to 



4 = -IGAtt - 4>1 



/i 



tan 77 + 



hi 



+ (c.c.) 



(A.7) 



(.1.8) 



(A9) 



where "(c.c)" means the complex conjugate quantity. Let us introduce the complex 
quantity 



P = l>l+ IV2 



(AlO) 



^/T^ V 1 + 

If 1? changes from— tt to tt, the changes of arguments of complex quantities tan | + 
and tan f — are respectively 

A arg(tan — + i^) = — 7rsgn(Imz/), A arg(tan — — i^) = 7rsgn(Imi/) (^-H) 
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Then 



Let 



-16A7r + 8 Ani 



/i 



= - (c.c.) 



sgni/2 



II — III + iii2 — cos(Ai + 1X2) — cos A 
where //i, //2, Ai, A2 are real quantities. Then according to (A. 10) 



= i^l + iV2 



' 1 — cos A 



tan 



A 



sin Ai 



/ 1 + cos A 
It follows from (A. 13) 

111 — cos Ai cosh A2 



sinh A2 



2 cos Ai + cosh A2 cos Ai + cosh A2 



II2 — sin Ai sinh A2 



Or 



sin^ Ai = |g 



1/2 



III + III -I 



cos^ Ai = 



sinh A2 
where 



III + III -I 



+ \q 



1/2 



cosh A2 ~ 



/^l + /^2 + 1 

2 

III + III + I 



1 1/2 



+ \q 



1/2 



f iil + iil + i 



{A.12) 
{A.13) 

{A.U) 
{A.15) 

(A16) 
(A17) 

{A.18) 



and the square root of the module is supposed to be non-negative. 

All functions (A. 16), (A. 17) are non-negative. They vanish only at the folowing 
conditions: sin^Ai = at //2 = OA | /Xi |> 1, cos^Ai = at /xi = 1, sinh^A2 = at 
;U2 = OA I |< 1. It means that sinhA2 is a single- valued function on the complex 
plane 11 with the cut A*2 = OA | /xi |< 1, and the sign of 1/2 



sgni/2 = sgn 



sinh A2 



cos Ai -I- cosh A2 



= sgn(sinhA2) 



(A19) 



is the same on the whole complex plane because | cos Ai |< cosh A2. cos Ai may 
coincide with — coshA2 only if sin^Ai = 0, and hence H2 — OA \ iii |> 1. In this 
case cosh^A2 = (/i? + l)/2+ | (/x^ - l)/2 |> 1 for iil > 1. cosh^ A2 = 1 only if 
/Li2 = A = ±1, but these points /x = 1 and /x = —1 lie on the cut and are 
supposed to be inaccessible. Thus 1^2 has the same sign on the whole complex plane 
IJ, with the cut 112 — OA 
represented as follows 



i - (c.c.) 



Ill \< 1. Expression in the brackets of (A.12) can be 

cosAi cosAi ,7r ^. .tt 

; =^ = tan( A — tan( A = 

sinAi smAi ^2 ^ ^2 ^ 



2 sinh A2 cosh Ao 



cosh A2 + sinh A2 — cos^ Ai -|- sin Ai 



{A.20) 
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Using (A. 16), (A. 17), one obtains 

cosh^ A2 + sinh^ A2 — cos^ Ai + sin^ Ai = 4 



1/2 



(A21) 



sinh^ A2 cosh^ A2 = ( d±A^\ /l±ifi 



(A. 18) is positive everywhere on the plane with the cut = OA | //i |< 1 and 



{A.22) 



h = -16A7r + 4>l7r 



-sgn(sinhA2) 

(^.23) 



is a single-valued function of arguments = ^ /^a and = r] /^a ■ 
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